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Clustering of luminous red galaxies II: 
small scale redshift space distortions 
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ABSTRACT 

This is the second paper of a series where we study the clustering of LRG galaxies 
in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galaxies 
covering over 1 Gpc^ /h^ for 0.15 < z < 0.47. Here we focus on modeling redshift space 
distortions in ^(ct, tt), the 2-point correlation in separate line-of-sight and perpendicular 
directions, at small scales and in the line-of-sight. We show that a simple Kaiser model 
for the anisotropic 2-point correlation function in redshift space, convolved with a 
distribution of random peculiar velocities with an exponential form, can describe well 
the correlation of LRG at all scales. We show that to describe with accuracy the so 
called "fingers-of-God" (FOG) elongations in the radial direction, it is necessary to 
model the scale dependence of both bias b and the pairwise rms peculiar velocity (J12 
with the distance. We show how both quantities can be inferred from the ^(ct, tt) data. 
From r ~ 10 Mpc/h to r ~ 1 Mpc/h, both the bias and au are shown to increase by 
a factor of two: from 6 = 2 to 5 = 4 and from 0-12 — 400 to 800 Km/s. The later is in 
good agreement, within a 5 percent accuracy in the recovered velocities, with direct 
velocity measurements in dark matter simulations with O^ = 0.25 and (78=0.85. 



1 INTRODUCTION 



The luminous red galaxies (LRGs) are selected by color and 
magnitude to obtain intrinsically red galaxies in Sloan Digi- 
tal Sky Survey (SDSS) (Eisenstein et al 2001). These galax- 
ies trace a big volume, around lGp(?h~'^ , which make them 
perfect to study large scale clustering or smaller scales with 
larger statistics. In this paper we focus on the 2-point cor- 
relation function ^(cr, tt) on the smaller scales, where the 
non-linear bias and the random peculiar velocities can com- 
plicate the analysis. Our motivation is to provide a model 
that can explain the observed ^(ct, tt) of the LRG galaxies. 
Such a model does not exist right now and we are not aware 
of any attempt to reproduce the ^{g, n) in the detail we will 
explore here. LRG galaxies seem to display larger "Fingers- 
of-God" (Jackson 1972) than regular galaxies, is this evi- 
dence for larger velocities? If so, is this evidence that LRG 
trace stronger gravitational potentials? We will show that a 
simple Kaiser model for the anisotropic 2-point correlation 
function in redshift space, convolved with a distribution of 
random peculiar velocities with an exponential form, can 
describe well the correlation of LRG at all the scales. To de- 
scribe with accuracy the so called "fingers-of-God" (FOG) 
elongations in the radial direction, it is necessary to model 



the scale dependence of both bias b and the pairwise rms 
peculiar velocity (T12 with the distance. 

In Paper I (Cabre and Gaztafiaga 2008) of this series, 
we have analyzed larger scales and present the basis for this 
work, including more detailed theory, error and systematic 
effects. 

On large scales, the density fluctuations are small enough 
to be linearized, and can be used to constrain cosmological 
parameters, since we can assume that the clustering is well 
described by (linearly biased) dark matter (see Paper I). On 
smaller scales, we can learn about the relation of galaxies 
to dark matter through the biased clustering of halos. This 
range of distances can be fitted by a power law, but there 
are small deviations that can be understood in the theory 
of the halo occupation distribution. The transition between 
galaxy pairs of the same halo and galaxy pairs that belong 
to different halo, occurs around IMpc/h. When moving to 
scales smaller than IMpc/h, in the 1-halo term, we can see 
processes more complex that modify the galaxy clustering, 
such as dynamical friction, tidal interactions, stellar feed- 
back, and other dissipative processes. 

We use the most recent spectroscopic SDSS data re- 
lease, DR6 (Adelman-McCarthy et al 2006), to perform 
the study of small scales in the anisotropic 2-point correla- 
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tion function, and its derivatives projected correlation func- 
tion and real-space correlation function. LRGs are supposed 
to be red old elliptical galaxies, which are usually passive 
galaxies, with relatively low star formation rate. They have 
steeper slopes in the correlation function than the rest of 
galaxies, since they are supposed to reside in the centers of 
big halos, inducing non-linear bias dependent on scale, for 
small scales. 

The same LRGs (but with reduced area) have been 
studied from different points of view. Zehavi et al (2005) 
study LRGs at intermediate scales (0.3 to 40Mpc/h), where 
they calculate the projected correlation function, the monopole 
and real-space correlation function to study mainly linear 
bias, non-linear bias and the differences between luminosi- 
ties. They find that there are differences from a power law 
for scales smaller than IMpc/h and they find no strong evo- 
lution for LRG. At smaller scales, Eisenstein et al (2005) did 
a cross- correlation between spectroscopic LRG with photo- 
metric main sample in order to reduce shot-noise in small 
scales clustering. They conclude that the clustering is higher 
for most luminous galaxies. Moreover, LRG have a scale de- 
pendent bias on luminosity, while for normal galaxies the 
luminosity bias is scale independent. LRG galaxies are sur- 
rounded by other red galaxies near them, which seem to be 
approaching the center LRG galaxy. Finally, Masjedi et al 
(2006) deal with very small scale clustering to scales smaller 
than 55" by cross- correlating the spectroscopic LRG sample 
and the targeted imaging sample and find that the correla- 
tion function from 0.01-8Mpc/h is really close to a power law 
with slope -2, but there are still some features that diverge 
from the power law. 

The small scale slope depends on the interplay between 
two factors which control how the correlation function of 
galaxies is related to that of the underlying matter : the 
number of galaxies within a dark matter halo (HOD) and 
the range of halo masses which contain more than one galaxy 
(Benson et al 2000). LRGs are to be found in halos with 
the median of the distribution occurring at 3 10 Mq/Zi, 
estimated using weak lensing measurements (Mandelbaum 
et al 2006). Almeida et al (2008) find with simulations that 
25% of LRGs at z = 0.24 are satellite galaxies, which play 
an important role in the form of the small scale correlation 
function, as well as the pairwise velocity dispersion, and also 
provide information about galaxy formation and evolution. 
Zheng et al (2008), who study HOD in LRG clustering, find 
that the satellite fraction is small (5-5% for Mg < —21.2) 
and decreases with LRG luminosity. 

Slosar et al (2006) show that pairwise velocity distribu- 
tion in real space is a complicated mixture of host-satellite, 
satellite-satellite and two-halo pairs. The peak value is reached 
at around IMpc/h and does not reflect the velocity disper- 
sion (712 of a typical halo hosting these galaxies, but is in- 
stead dominated by the sat-sat pairs in high-mass clusters. 
Different groups have found this dependence of ai2 on the 



scale (Zehavi et al 2002, Hawkins et al 2003, Jing and Borner 
2004, Li et al 2006, Van den Bosch 2007, Li et al 2007). They 
also find that cri2 depends on luminosity, with higher values 
for both large and small luminosities, but this tendency can 
not be predicted by the halo model. Li et al (2006, 2007) add 
that redder galaxies have stronger clustering and larger ve- 
locities at all scales, so they move in strongest gravitational 
fields. Tinker et al (2007) use the halo occupation distribu- 
tion framework to make robust predictions of the pairwise 
velocity dispersion (PVD). They assume that central galax- 
ies move with the center of mass of the host halo and satel- 
lite galaxies move as dark matter. The pairs that involve 
central galaxies have a lower dispersion, so the fraction of 
satellites strongly influences both the luminosity and scale 
dependence of the PVD in their predictions. At r between 
1 and 2 Mpc/h, the PVD rapidly increases with smaller 
separation as satellite-satellite one-halo pairs from massive 
halos dominate. At r < 1 Mpc/h, the PVD decreases to- 
wards smaller scales because central-satellite one-halo pairs 
become more common. At r > 3Mpc/h, the PVD is donn- 
nated by two-halo central galaxy pairs, and reach a constant 
value. 

LRGs have also been analyzed at higher redshifts (z=0.55) 
with the 2dF-SDSS LRG and QSO Survey (2SLAQ, Cannon 
et al 2007). Ross et al (2007), Da Angela et al (2008) and 
Wake et al (2008) analyze the redshift distortions in the 
LRGs and quasars for this catalog. For part of our analy- 
sis we have followed the method explained in Hawkins et al 
(2003), an extensive analysis of the redshift distortions in 
the 2dF catalog. 

Here we define the parameters that we assume during all 
this work, which are motivated by recent results of WMAP, 
SNIa and previous LSS analysis: Us — 0.98, fit = 0.045, 
h — 0.72 and fiat geometry. Unless otherwise said, we use 
^m = 0.25 and as — 0.85, both obtained from paper I. We 
will use the power spectrum analytical form for dark matter 
by Eisenstein and Hu (1998), and the non- linear fit to halo 
theory by Smith et al (2003). 



2 MODELING REDSHIFT-SPACE 

We follow the modeling given in more detailed in Paper 
I (Cabre and Gaztaiiaga 2008). Here we just show the 
main equations related to this work. In the large-scale lin- 
ear regime and in the plane-parallel approximation (where 
galaxies are taken to be sufficiently far away from the ob- 
server that the displacements induced by peculiar velocities 
are effectively parallel), the distortion caused by coherent 
infall velocities takes a particularly simple form in Fourier 
space (Kaiser 1987); 



Ps{k) = (l + pdfP{k). (1) 
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where P{k) is the power spectrum of density fluctuations S, 
H is the cosine of the angle between k and the hne-of-sight, 
the subscript s indicates redshift space, and /3 is the growth 
rate of growing modes in linear theory. 
The observed value of /3 is 

/(n„0 _ IdlnD 



/3 = 



(2) 



b h din a 

where b is the bias between galaxies and dark matter, /(fi,„) 
and D(yirn) are the linear growth density and velocity fac- 
tors. Here we use Hamilton (1992) who translated Kaiser 
results into real space ^'(cr, tt) (see Eq.(8) in Paper I). We 
then convolve it with the distribution function of random 
pairwise velocities, /(«), to give the final model ^(cr, tt) (Pee- 
bles 1980): 

/OO 
i'{o,T,^v/H{z)/a{z))f{v)dv (3) 

-etc 

We represent the random motions by an exponential 
form (Szalay et al, 1998), 



/(«) 



1 ( V2\v\ 

: exp 



(4) 



where 012 is the pairwise peculiar velocity dispersion. 

It is well known that 0-12 depends on the real sepa- 
ration between particles r = \/ cr'^ + tt'^^^i, where iVreai = 
Tr~v/H{z)/a{z). The pairwise velocity dispersion is roughly 
constant for r > 5Mpc/h, but increases towards smaller val- 
ues of r. To use Eq.(13| and Eq.Q with (T12 depending on 
r, we do the following. For each value of the velocity v in 
the integral, we estimate the real distance r = \/o-^ + '""reai 
and the value of (T12 = (Ti2(r) that enters in f{v). We have 
checked that for a > 2Mpc/h, this exact method is very sim- 
ilar to assuming that a\2 is fixed for each a and it is constant 
along the line-of-sight tt, so we can change cri2 depending on 
(J, rather than r. This has some practical advantages over 
the exact method. 

We define the multipoles of ^((j, tt) as 

Us) = ^^l^\{<^,n)Pdfi)df,. (5) 

where fi is cosine of the angle to the line-of-sight tt. The 
monopole ^(s) = Co(s) is also called the redshift space cor- 
relation function. The real-space correlation function can be 
estimated from the projected correlation function, H(cr), by 
integrating the redshift distorted ^(a, vr) along the line-of- 
sight n: 



H(a) = 2/ C(a,^)d7r 



(6) 



Davis and Peebles (1983) show that E(cr) is directly 
related to the real-space correlation function. 



S(a) = 2 r 



r^{r)dr 
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(7) 



It is possible to estimate ^(r) by directly inverting ^{(t) 
(Saunders et al 1992). We can write Eq.lJT]) as. 



cw = -- 



(d=.(G)/da) 



do. 



(8) 



''^ Jr (p'^ — r^) 2 

Assuming a step function for H(cr) = E.i in bins centered on 
(Ti, and interpolating between values. 



e(^o 






71" ^ CTj + l - CTj 



a, + J a 



(9) 



for r — ai. 

Once we recover the real-space correlation function, we 
can also estimate the ratio of the redshift-space correlation 
function, ^(s), to the real-space correlation function, ^(r), 
which gives an estimate of the redshift distortion parameter, 
/3, on large scales: 



f--f-f 



3 STUDY OF THE ERRORS 

A detailed account of errors is given in Paper I. Basically 
we use mock catalogs to estimate what we call the Monte 
Carlo (MC) errors. Mock catalogs are build out of very large 
numerical simulations run in the super computer Mare Nos- 
trum in Barcelona by MICE consortium (www.ice.cat/mice). 
The simulation contains 2048^ dark matter particles, in a 
cube of side 7680Mpc/h, Qm = 0.25, fit = 0.044, erg = 
0.8, Us = 0.95 and h = 0.7. We use both dark matter 
and groups. There is no bias in the dark matter mocks, so 



zf + 1 



n„(2)" 



= 0.62 (where n„(z) = f7™(l + z)7(r2„(l + 
Qm) and b — 1). In this paper, we use group mocks 
with M > 2.2 X 1O"M0 (b = 1.9, (3 = 0.25) at z = 0, which 
are similar to real LRG galaxies in its clustering properties. 

For the jackknife (JK) error, we obtain the different 
realizations directly from the data, dividing the catalog in 
M zones, and we consider that each realization is all the 
catalog except from one of these JK zones. In this case, as 
the realizations are clearly not independent, we multiply the 
covariance matrix by a factor {M — 1) to account for this 
effect (see Paper I for more details). 

To analyze data we prefer to use a model independent 
error, such as the JK error above. The MC error based on 
simulations depends on the model that we have used and in 
particular in the overall normalization, which is similar to 
LRG for group simulations but can have slight differences. 
We use the MC errors to probe the accuracy of JK errors in 
the same simulations. 

In Paper I (Cabre and Gaztaiiaga 2008), we present 
the errors in the redshift space correlation function ^{a, n) 
(bin=5Mpc/h), the monopole ^(s) and the quadrupole Q{s). 
Here we study the errors in the ^{a, n) (bin=lMpc/h and 
0.2Mpc/h), perpendicular projected function H(a) and in 
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Figure 2. We compare different error estimators for ^(cr, n) with 
IMpc/h of bin for MICE group mocks with b = 1.9. We fix the 
perpendicular distance a and move along the line-of-sight tt as 
indicated in the figure. The errors are: MC (solid line), theory 
with (dashed) and without (dotted) the signal part of the error, 
and JK (big dots with errorbars). As we move to higher a (lower 
amplitudes), JK starts to fail, but the analytical and MC errors 
there agree well. 



JK error 



Analytical error 
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Figure 3. We compare different error estimators for ^(cr, tt) with 
0.2Mpc/h of binning for LRG real data: JK and analytical. As in 
the case of bin=lMpc/h, JK seems to work for a < 20Mpc/h. 



the real-space correlation function ^(r). We also use the 
simulations to test the validity of the 2-point correlation 
function model. 



3.1 Errors in the redshift space correlation ^(a, vr) 

We calculate the 2-point correlation function in redshift space 
^((T, tt) for each group mock, using a bin of IMpc/h. We 
can obtain a JK error for each mock, so a mean and dis- 
persion for all the mocks. We compare it to the MC error, 
and to the analytical form, which is described in detail in 
paper I. We propose the error to have the following form 
A^ = A^shot-noise+A^signai, with two arbitrary coefhcients 
Onoise and Osignai , SO that: 



^S — C^noise ^i,Poisson ~r Oisignal S V j 

The comparison for the three kind of errors can be seen in 
Fig[TJ where it is also plotted the JK error for real LRG 
galaxies, and the corresponding analytical error. The ampli- 
tude and shape of the error is similar in the group mocks 
and in LRG data, so we assume that our conclusions de- 
rived from the simulations can be translated to real LRG 
data. In order to see the differences in the mock errors with 
more detail, we fix the perpendicular distance a for different 
cases and move along the line-of-sight tt in Fig[2] The JK 
error starts to deviate slightly from MC for a higher than 
20Mpc/h, where the shot noise analytical form works per- 
fectly. For small a and tt, the signal part of the analytical 
error helps to fit to MC error, but it seems that at these 
scales the best option is to use JK error. The covariance is 
lower than 0.2 for all the points, so it is nearly diagonal. In 
FigOwe compare the JK error and analytical error for real 
data for a bin of 0.2Mpc/h and we find similar conclusions. 
We have done the same analysis with dark matter simula- 
tions and they also work well, as expected. Note that con- 
trary to what we found in paper I for 5Mpc/h binning, for 
smaller bins the shot noise model matches the expectations 
with a„oise ~ 1, rather than Onoise = 1.4. We believe that 
this could be related to the smaller covariance in 0.2Mpc/h 
and IMpc/h binning. 

3.2 Errors in the projected correlation H(cr) 

We calculate the projected correlation function H((7) for each 
group mock integrating through tt the redshift space corre- 
lation function ^(cr, tt), and we also calculate the JK error. 
Then we look at the difference between the MC dispersion 
and the mean over all the JK errors. In Fig|4l top panel, we 
see the mean H((t) over all the simulations (solid line with 
errors), and over-plotted in blue big dots the H((t) for the 
real LRG data. In the bottom panel we show a comparison 
of the different errors in the simulations (MC dotted line, 
and JK solid line with errors) and the error JK in the LRG 
data (big dots). Errors in MC and JK in simulations coincide 
at the scales where we can use the S,{a) (below 30Mpc/h). 
The JK error in LRG is also similar to the errors in the sim- 
ulations. We have done the same analysis with dark matter 
simulations and they also work well, as expected. We con- 
clude that using JK errors gives a good approximation to 
the true MC errors. 



3.3 Errors in the real-space correlation ^(r) 

We now calculate the real-space correlation function ^(r) 
from the projected correlation function H(a) (Eq.(|9])). In 
FigE top panel, we see the correlation function obtained 
from deprojecting H(cr) with errors (solid line) compared to 
the real-space correlation function which we can calculate 
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Figure 1. We compare different error estimators for g(o-, tt) with binning IMpc/h for MICE group mocks with b = 1.9: Monte Carlo, 
mean of all jackknife and analytical form. The last two plots refer to real LRG data: JK error and analytical form. Errors work well in 
the simulations, and they are similar in amplitude and shape to LRG real errors. 
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obtained in LRG data, with a similar bias than the group 
simulations, as expected {b was found to be fe ~ 2 in Paper 
I). On the bottom panel of FigO the plot compares errors 
in MC (solid line with errors) and JK (dotted line) case, and 
they are very similar at small scales. We have over-plotted 
the error JK in LRG data (big dots). We have also done a 
similar analysis with dark matter mocks and they also show 
a good agreement between JK and MC errors. As in the 
previous case, we conclude that using JK errors give a good 
approximation to the true MC errors. 



CT(Mpc/h) 



10.0 




CT(Mpc/h) 



Figure 4. Top: S(ct) for the group simulations calculated from 
^(c, 7r) (solid line with errors) and the value for LRG data (big 
dots). Bottom: Errors in S((t) for MC simulations (dotted line), 
JK simulations (solid line with errors) and JK for real LRG data 
(big dots) 



using the simulations in real space, without redshift dis- 
tortions (dotted line). We can see in Fig[5]how we recover 
well 4(r) for small scales (below 40Mpc/h). Over-plotted in 
blue (big dots), we see the real-space correlation function 



3.4 Validity of the models 

Besides studying errors, we also use the simulations to test 
the methods that we will apply to real data (LRG). In top 
panel of Fig[5]we showed that we can recover the real-space 
correlation function. Now we want to see if we can model the 
redshift-space correlation function with the simple model ex- 
plained in !j2l when taking the approximation of ai2 constant 
along the LOS and varying it for a lower than 5Mpc/h. In 
Fig[6] we plot ^{a, tt) at small scales for the mean of the 
MICE dark matter mocks. We over-plot in solid lines our 
model with a varying pairwise velocity dispersion. We based 
this model on the pairwise velocity dispersion ai2{r) esti- 
mated from the velocities in the simulations. Fig[7]shows the 
measured pairwise velocity dispersion cri2(r). In our model, 
for a larger than 5Mpc/h we use the effective ai2 = 400fcm/s 
(asymptotic value at large scales). For a lower than 5Mpc/h 
we use a different cri2 as given by FigO With this simple 
approximation, we reproduce very well the observed FOG 
as shown in Fig|6] We have done this study with dark mat- 
ter mocks since we know the variation of (Ji2 exactly, while 
when we construct the group mocks, we loose the pairwise 
velocity at small scales due to the way we identify groups. 
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Figure 5. Top: ^(r) for the group simulations calculated from 
H(r) (solid line with errors) and the value for LRG data (big 
dots). Bottom: Errors in ^(r) for MC simulations (dotted line), 
JK simulations (solid line with errors) and JK for real LRG data 
(big dots) 



4 THE DATA: LUMINOUS RED GALAXIES 

SDSS luminous red galaxies (LRGs) are selected on the basis 
of color and magnitude to have a sample of luminous intrin- 
sically red galaxies that extends fainter and farther than the 
SDSS main galaxy sample. Eisenstein et al (2001) give an 
accurate description of the sample. In Paper I we give more 
details of our selection. Here we just give a summary. 

We k-correct the r magnitude using the Blanton pro- 
gram 'kcorrect' q We need to k-correct the magnitudes in 
order to obtain the absolute magnitudes and eliminate the 
brightest and dimmest galaxies. We have seen that the pre- 
vious cuts limit the intrinsic luminosity to a range —23.2 < 
Mr < —21.2, and we only eliminate from the catalog some 
few galaxies that lay out of the limits. Once we have elimi- 
nated these extreme galaxies, we still do not have a volume 
limited for high redshift galaxies, but we suppose that the 



http://cosmo.nyu.edu/blanton/kcorrect/kcorrect_help.html 
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Figure 6. ^(cr, w) for the MICE dark matter simulation (as col- 
ors). The contour colors are -0.05, -0.01, -0.005, -0.001, 0, and 
0.001 to 20 with 20 equally spaced logarithmic bins. In solid 
lines we have plotted the model for the input parameters of 
the simulation. For a larger than 5Mpc/h we use the effective 
(TJ2 = 400A;ni/s (asymptotic value at large scales) to model ^(cr, it) 
while for a lower than 5Mpc/h we use a different cri2 as seen in 
Fig[7l for each cr constant along the LOS n 
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Figure 7. Dispersion in the pairwise velocity distribution (T12 
estimated from the velocity field in the MICE dark matter simu- 
lations as a function of the distance between particles. 



variations in luminosity just change the overall shape in the 
clustering. 

We have masked the catalog using at the first step the 
photometric DR6 mask, based on the number of galaxies 
per pixel. In previous works we saw that the mask that we 
obtain statistically by dropping out the pixels with small 
number of galaxies gives the same correlation function that 
the one obtained by extracting the polygons masked by the 
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SDSS team. After that, we compare our masked catalog to 
the LRG spectroscopic catalog, and we extract the galaxies 
that lay outside from "good" plates. 

This rough mask could imprint spurious effects at very 
small scales, but we are not interested in these scales where 
fiber collisions in the redshift catalog are limiting our analy- 
sis, for distances less than 55arc sec, less than 0.3Mpc/h at 
the mean redshift of LRG data, z=0.35. We obtain 75,000 
galaxies for the final catalog, from z=0.15 to z=0.47. The 
area of the data used is around 13.5% of the sky. See Fig.2 
in Cabre and Gaztanaga (2008) for a plot of the mask. Re- 
cently we have used another mask, provided by Swanson et 
al. (2008), which is in a readily usable form, translating the 
original mask files extracted from the NYU Value-Added 
Galaxy Catalog (Blanton et al. 2005), from MANGLE into 
Healpix format (Gorski et al. 1999), and results are very 
similar. We have also done the analysis with the recent data 
release DR7, the final one in SDSS, and results are very 
similar at small scales. 

We use the ^ estimator of Landy and Szalay (1993), 



^(a,^) = 



DD - 2DR + RR 
RR 



(12) 



to estimate the 2-point correlation function in redshift 
space, with a random catalog 20 times denser than the SDSS 
catalog. The random catalog has the same redshift distribu- 
tion as the data, but smoothed to avoid the elimination of 
intrinsic correlations in the data, and also the same mask. 
We count the pairs in bins of separation along the line-of- 
sight (LOS), n, and across the sky, a. The LOS distance is 
just the difference between the comoving distances in the 
pair. The perpendicular distance is cr = \/s'^ — tv^, which 
corresponds approximately to the mean redshift. 



5 ANALYSIS 

As we explore small scales, we will show in this section 
that we encounter the following problems. First of all, the 
bias becomes clearly dependent on the scale for distances 
smaller than about 10 Mpc/h, because LRG are galaxies 
highly biased so they only keep the linear bias constant at 
large scales. Secondly, the model that we are using assumes 
that the pairwise velocity dispersion is independent of scale, 
which is not a good approximation for small scales as we 
have seen in the simulations. We can arrange this by using 
different ai2 for each real distance r, as explained in Section 
[2l Here we want to check if we can infer the correct varia- 
tion of (Ti2(r-) with scale using the ^(ct, vr) data. As a test, 
we will compare the inferred values of cti2 (r) with the direct 
measurement from the velocity field in the simulation. 
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Figure 8. l;{r) (blue dots) and ^(s) (orange dots). We see clearly 
how the redshift-space correlation function is the real-space corre- 
lation function biased by a constant factor that represents gravi- 
tational infall (dependent on /3 in Kaiser approximation) at large 
scales above 4Mpc/h (although ^(r) is overestimated for scales 
larger than 30Mpc/h due to the precision in the calculation). 
However, for small scales, the redshift-space ^(s) is strongly sup- 
pressed compared to ^(r) due to random peculiar velocities. 



5.1 Real space correlation 

In Fig|5]we show the resulting real-space correlation function 
which we have calculated using Eq.® (in blue) and over- 
plotted the monopole in redshift space (in orange). At inter- 
mediate scales, from 5 to 30 Mpc/h (the top value limited by 
the method to obtain ^(r)), f (s) is equal to ^(r) but biased 
by a constant factor, as in Eq. (|10|l . due to Kaiser effect. We 
had shown this constant bias at intermediate scales (around 
lOMpc/h) in Fig. 8 of paper I (Cabre and Gaztafiaga 2008), 
where we plot the ratio of the correlation function in red- 
shift space and real space iis)/^{r). At large scales, we can 
cissociate it to a function of f3 obtained in Paper I at large 
scales. The agreement is excellent, which provides a good 
consistency check for our results. The difference between the 
real and redshift space correlation function at small scales 
is primordially due to the random peculiar velocities. 

5.2 Power law fit 

We next fit our estimation of the real space correlation func- 
tion to a power law, ^(r) = {r/ro)~''°, from IMpc/h to 
15Mpc/h. At scales smaller than IMpc/h the fit is not good, 
it is no longer a power law. In Fig [9] we show a fit for ro and 
7o. In Fig llOl we show the measured real space correlation 
function, the best fit power law model (red) and the dark 
matter non-linear correlation function obtained from best 
parameters in Paper I which works well for a linear bias, 
ie on large scales. The large scale and small scale fittings 
models agree well at large scales, as can be seen in the plot, 
and the correlation function does not follow a power law for 
distances smaller than ^ IMpc/h, where we can probably 
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Figure 9. Best fit of 70 and tq of ^(r) to a power law model 
5(r) = (r/ro)"™, l-15Mpc/h. Dotted and continuous contours 
correspond to 1 ad 2 degrees of freedom with 1,2 and 3-sigma 
confidence. 
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Figure 10. Observed ^(r) (symbols with errorbars), best fit to 
the power law ^(r) = (r/ro) '" (red), and dashed over-plotted 
the best real space correlation function for large scales, assuming 
a constant bias 



see the transition from the one-halo to the two-halo term 
(Scoccimarro 2004). 

These results are similar to other studies. Zehavi et al 
(2005) did a similar analysis with a previous SDSS spectro- 
scopic data release (35000 LRGs) at intermediate scales from 
0.3 to 40Mpc/h. We have doubled the number of LRGs and 
our results agree with them for the monopole, the projected 
correlation function, and the obtained real-space correlation 
function, with the same main conclusions. Also Eisenstein et 
al (2005) , in a study of small scales (0.2-7Mpc/h) using the 
cross-correlation between spectroscopic LRG with the main 
photometric sample, remark that ^(r) can not be explained 
with a power-law fitting. However, Masjedi et al (2006) have 
obtained the correlation function at very small scales (0.01- 
8Mpc/h) and have found that, although with some features 
diverging from a power law, all the range is really close to a 
5(r) oc r^^ . 

5.3 Non-linear bias 

We have checked in our simulations with halos that the non- 
linear bias typically follows a power-law (at least for dis- 
tances larger than 0.3Mpc/h), which has a different slope 7;, 
depending on the halo mass. In general this could depend 
on other parameters concerning galaxy formation. LRG are 
assumed to be red galaxies that trace halos of 10^^ M©, but 
there is a wide range of halos masses, and the non-linear 
bias shows us these properties. Here we estimate the bias as 



h{r) = ^i{r)/i{r)DM = bhr.i{r) (13) 

where we have separated the bias &(r) into a constant 
scale independent term, b, on large scales and a function of 
scale b„i{r) which goes to unity on large scales. When b„i{r) 
reaches unity, we assume that the bias is linear from there 



on. We define a parameter ri, for the power law 6„i(r) = 
(j'/ri,)~'"' which shows the scale from which non-linearity 
begins to be important. The value of rj should coincide ap- 
proximately with the correlation length, where the real space 
correlation function is unity. 

To compute the non-linear bias b(r) from the above 
equation, we need the dark matter correlation function and 
the linear bias &, which we have taken from the fitting at 
large scales done in Paper I. We have calculated the bias 
for all the values of 57^ and different amplitudes. Then we 
marginalize over them. In Fig llll we see the contours for 
ri, and 7t, (top panel), and the best fit (red in the bottom 
panel). This fit to the bias can explain the differences seen 
previously between the correlation function and a power 
law for scales smaller than IMpc/h. At scales smaller than 
0.3Mpc/h, the real space correlation function turns down 
due to fiber collisions (Masjedi et al 2006). In detail, we 
see in Fig llll a feature in the bias between 1 and 2 Mpc/h, 
indicating that LRGs galaxy bias is not completely smooth. 
We think that this feature is due to the range of halo sizes 
of our LRGs, which makes it difficult to predict exactly the 
transition point from the 1-halo to the 2-halo term (Scoc- 
cimarro 2004). If galaxies are residing within dark matter 
halos then the clustering of the galaxies on scales larger than 
halos is determined by the clustering of the dark matter ha- 
los that host them, plus statistics of the occupation of halos 
by galaxies. For larger scales than ~ 2Mpc/h (the biggest 
halos), the clustering comes entirely from LRGs that reside 
in different halos, while for smaller scales, the clustering can 
come from galaxies in different halos or galaxies in the same 
halo until it is reached a minimum size of halos (see Masjedi 
et al (2008) for a more detailed explanation). 
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Figure 11. Top panel: Best fit to no linear bias b„i{r) (defined 
in the text) with a power law b„i{r) = (r/ri,) ">. Bottom panel: 
Non linear bias b^iir) (solid line with errors in gray) and best 
power law fit (red). We have also over-plotted in dashed line the 
bias obtained if we suppose that the galaxy correlation function 
is a power law. 
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Figure 12. Monopole ^(s) (dots) with errors (gray), best model 
assuming a constant ai2 (solid red) and model assuming variation 
in (Ti2 with errors (dashed blue), as in Fig llSI Bottom panel: 
Quadrupole ^2{s) (dots) as in the top panel. 



measurements. These differences indicate that ai2 is higher 
at smaller scales as shown in Fig[7]for simulations. 



5.4 Monopole and Quadrupole 

Once we have obtained the real space correlation function, 
we can look at the monopole ^(s), the quadrupole (,2{s) and 
also at ^{a, tt) in order to check the result. We see in Fie: ll2l 
the monopole ^(s) (top panel) and quadrupole ^2(s) (bot- 
tom panel) binning the distance with 0.2Mpc/h, and over- 
plotted in red solid line the theoretical model, which we 
have found integrating ^{a,n) and assuming a constant ai2 
(derived using the normalized quadrupole Q(s) in Paper I). 
The prediction uses the model explained above and in Pa- 
per I, where the shape of ^{a, vr) is given by the real-space 
correlation function and the parameters describing the ve- 
locity distortions. The monopole directly measured in the 
data (dots with shaded region in top panel of Fig[T2} does 
not agree with the model (red line), which is lower for scales 
smaller than 3Mpc/h. The same happens to the quadrupole 
(bottom panel), where the model is also higher than the 



5.5 Recovering ai2 from simulations 

We now recover the dispersion of pairwise velocities ai2 at 
each perpendicular distance a assuming that cti2 remains 
constant along the LOS tt. In reality, ai2 is a function of the 
real distance r, ai2{r). At each n, a, the correlation is a con- 
volution of all the real distances. We could in principle try 
to fit all the tt-ct plane for different shapes of (Ti2(r), but this 
involves too much freedom and such a fit is not feasible in 
practice. The alternatives are to do the fit in Fourier space, 
where there is no dependence of ai2{k) along the LOS, or 
to use the approximation ai2{r) ~ (Ti2((t) in the zone of 
the n — a plane where this is a good approximation. We 
have studied the differences between ^{a, tt) obtained from 
either ai2{r) or from cri2 (cr) to explore this later possibility. 
The two models differ, more or less strongly depending on 
the case, at small a and large n. For our dark matter sim- 



© 0000 RAS, MNRAS 000, 000-000 



10 Cabre & Gaztanaga 



1200 


- 




1000 


- 




^ 800 




. • 


6 

^ 600 


- 




*= 400 


- 




200 


- 













CT(Mpc/h) 

Figure 13. We plot the mean recovered (712 (r) {iTmax = 
5Mpc/h) for a bin 0.2Mpc/h (black circles with gray zone as er- 
rors) and for IMpc/h (red circles with dashed zone errors), com- 
pared to the original ai2{r) (dashed line) calculated in the real 
space simulation. We take each mock redshift space ^(cr, 7r) and 
model it following Eq. ((Sjl, from where we obtain the dependence 
of (TJ2 on the distance. 



ulations, the differences are smaller than in the LRG case, 
where they can be large enough to bias the final result. The 
best option, then, is to fit for each perpendicular distance 
o"i2(o") using ^((7, tt) up to a maximum nmax = 5—lOMpc/h, 
where both models agree well. If we go further in tt, the re- 
sult is biased to lower values of ai2 . 

In order to test the method to recover ai2{r) we have 
used our dark matter mocks, where we trust velocities to be 
accurate. We take each mock and recover CT12 (r) for nmax = 
5, 10, 18, 30Mpc/h. We have used a bin of 0.2Mpc/h for a 
= 0.3 - 5Mpc/h and a bin of IMpc/h for cr = 1 - lOMpc/h. 
In Fig llSI we plot the mean recovered (Ti2(cr) for 0.2Mpc/h 
(black circles with gray zone as errors) and for IMpc/h (red 
circles with dashed zone errors), compared to the original 
(712 (r) calculated directly from velocities in the simulation 
(dashed line). We use nmax = bMpc/h for this plot, but 
we obtain very similar results when we go further in n. 
The main conclusion of this analysis is that we recover cor- 
rectly the random dispersion of pairwise velocities, and sec- 
ondly that the method to derive cri2 works at least until 
nmax ~ 30Mpc/h, for dark matter. This method is good to 
obtain the variation of cri2 at small scales. At larger scales 
(around 5Mpc/h), ai2 reaches a constant, and we see that 
in our study the estimated (T12 is slightly biased to higher 
values. At these distances, variations on cri2 do not change 
much the model for ^{a,n), so this is not a problem for 
comparison to data. We should take into account that we 
are using a simplified model which assumes that pairwise 
velocities are exponential. This is a good approximation for 
small scales, as we see from the simulations, but the real 
dispersion is in fact skewed and it is not perfectly exponen- 
tial, so we can see slight variations from the real dispersion. 



However, to recover an unbiased result for ai2 from ^((j, tt), 
the most important point is to take into account the infall 
velocities in the model (Scoccimarro 2004). Here we include 
infall velocities through the parameter 13. 

We next calculate the x^ for each mock in the cen- 
tral region of the plane tt — a (defined buy the contours 
of ^(o", tt) > 1.0,0.57,0.43,0.24) and we obtain a mean and 
dispersion of the reduced x^ (ie the x^ divided by the num- 
ber of degrees of freedom). The distribution of reduced x^ 
values from the different mocks can be seen in Fig |14l the 
mean is equal 1 and this means that the model works per- 
fectly well for our simulations of dark matter. But note how 
the dispersion in the distribution of x^ is quite broad, there 
is therefore room for some divergences away from 1. 



5.6 Measurement of ai2 as a function of scale 

Once we have showed that we can recover ai2{r) and that 
this simple model works for dark matter, let's do the same 
for the LRG galaxies. In order to calculate the model ^(a, n) 
at each a and vr, we need to assume the distortion param- 
eter /3 and the real space correlation function £,{r). We use 
j3 = 0.34, found in Paper I of this series. We also use the 
real space correlation function ^(r) from the integration of 
^((T, tt) along the line-of-sight, see Eq.(|8} and Fig[8] As men- 
tioned above, the real space correlation ^(r) is not well recov- 
ered above 30Mpc/h, where we will use instead the measured 
monopole corrected by the distortion bias factor, Ea. (|10p . 
At large scales, the real space and the redshift space corre- 
lation function are almost linearly biased, except from the 
BAO peak where there are some smaller non-linear effects 
again, but these do not have an effect in the modelization of 
5((T, tt) at small scales, so this is a good approximation. We 
calculate 0-12(0-) up to n^ax ~ 5, 10, 15, 20, 25Mpc//i. For 
values above nmax =10Mpc/h, we find that 0-12(0) starts to 
be biased low for small values of o with respect to the true 
values of oi2(r). For higher values of o, we can fit until a 
larger nm.ax, reducing the errors. We use the same binning 
as in the simulations and the resulting 012 (r) is plotted in 
Fig |15l where we compare it to the true value from velocities 
in simulations (solid line). Both results have approximately 
the same amplitude, that depends strongly on the cosmolog- 
ical parameters. If there is no velocity bias, then our universe 
must be similar to our simulations. 

Now we can recalculate the monopole and quadrupole 
including the variation of 012 with scale. We obtain a good 
fit to the results, as we can see in Fig |12l ( as a blue dashed 
line) compared to the same prediction assuming a constant 
012 (solid red line). 

As mentioned in !jl] Slosar et al (2006) show that pair- 
wise velocity distribution in real space is a complicated mix- 
ture of host-satellite, satellite-satellite and two-halo pairs. 
The peak value is reached at around IMpc/h and does not 
reflect the velocity dispersion of a typical halo hosting these 
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Figure 14. Distribution of x^ obtained from the mocks when we fit each ^{a, tt) with the dark matter model which includes the correct 
cri2(r). Different panels are showing the zone where we fit, from smaller to larger, for ^{cr,iT) > 1., 0.57, 0.43, 0.24 
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Figure 15. cri2 vs distance for LRG galaxies (circles with error- 
bars), obtained from modeling ^(cr, tt) with Eq.Q. We compare 
it to the true value from velocities in the dark matter simulation 
(solid line). 



galaxies, but is instead dominated by the sat-sat pairs in 
high-mass clusters. Tinker et al (2007) uses the HOD model 
to explain that at r ^ 1 - 3 Mpc/h, the PVD rapidly in- 
creases towards smaller scales as satellite-satellite pairs from 
massive halos dominate. At r < 1 Mpc/h, the pairwise ve- 
locities dispersion decreases with smaller separation because 
central-satellite pairs become more common, but we do not 
see this tendency because we do not study such small scales 
due to fiber collisions. At r > 3 Mpc/h, the (T12 is dom- 
inated by two-halo central galaxy pairs. Their predictions 
agree well with our results in shape and amplitude (compare 
Fig llSl with their Fig.4 more luminous prediction). Note that 
the amplitude of the predictions can be increased or reduced 
easily by changing the cosmological parameters. 

Tinker et al (2007) also distinguish between ai2(r) and 
o"i2(o-). They found that ai2{r) is larger than cri2((T) (com- 
pare their Fig.4 with Fig. 6). We have shown with our sim- 
ulations that (T12 (a) can describe well (T12 (r) provided nmax 
is small enough. Different groups have found this depen- 
dence of (712 on the scale (Zehavi et al 2002, Hawkins et al 
2003, Jing and Borner 2004, Li et al 2006, Van den Bosch 
2007, Li et al 2007 and citations there). It is difficult to 
provide a detailed comparison with these previous results. 



This is because very different assumptions are used: the in- 
fall model, the modeling of the real-space correlation, the 
values of TTmax Or the use of Fourier versus configuration 
space analysis. However, these previous results seem to find 
lower values of ai2{cr) than our results, obtaining a max- 
imum pairwise velocities of around 600 km/s, rather than 
then 800 km/s that we find. We believe that this differ- 
ence is mainly caused by the methodology, which can make 
o"i2(o") lower that ai2{r) when large values of -Kmax are used. 
In our case we find with our modeling that smaller values 
of TTmax recover better the right values of ai2 (r) when the 
galaxy bias b is large, as in the case of LRG. These larger 
values at ai2{r) agree well with the dark matter prediction 
for Qrn = 0.25 and erg = 0.85, as found in paper I, so there 
is no need to postulate that LRG velocities are in any way 
different. 



5.7 Consistency of a vr model: FOG 

Now we look directly at ^(cr, vr) at small scales once we 
have all the parameters, to see if the model works when 
we separate vr and a, for all the angles, rather than just 
the monopole or quadrupole. We have used a binning of 
0.2Mpc/h for these plots, in order to see clearly the fingers 
of God, which are concentrated at very small a. First, we can 
see the detailed plot of the measurements ^(ct, tt) in Fig |16l 

In the next three figures 1171 [181 and 1191 we show the 
differences between the data and the model in three cases. 
Top panels from left to right show: the data as colors; the 
data and the model over-plotted as solid line; and the model 
as colors. Bottom panels show: the same as in the top but 
we have zoomed the a direction to see clearly the fingers of 
God. 

In Fig |17l we compare the data with a model that as- 
sumes linear bias b (found fitting large scales only) and a 
constant pairwise velocity dispersion of ai2 (obtained from 
the quadrupole Q(s)). We see clearly that this does not work 
well. The FOG (along the n direction) are too small compare 
to data and the correlation in the a direction has a different 
slope. This is partially due to the fact that the bias in the 
real data becomes non-linear on smaller scales. Part of the 
apparent lack of fingers of God in the models are corrected 
just adding the non-linear bias in the model. 
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Figure 17. Redshift space 2-point correlation function i;{cr, it) (for the main slice in redshift z=0. 15-0.47) modeled with a linear bias in 
the real-space correlation function, that we need to obtain ^(cr, tt) with Eq. JSj and an effective o-i2=380km/s, as find for LRG at large 
scales. Topidata (as colors), data (as colors) -|- model (solid line), model (as colors). Bottom: the same zoomed in a 
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which we correct with the third model in Fig |19l Here we 
include the variation of (Ji2 with scale, assuming that ai2 
is constant for large scales and changes for small scales, as 
found in Fig llSI Here we use the exact modeling of ai2 (r) as 
a function of r as explained in section [5] In order to see more 
clearly the validity of the model, we have plotted in Fig |20l 
the correlation function ^{a, n) along the line-of-sight n for 
different fixed a as indicated in the figure. Similar results 
are found for all other a values, which is studied in bins of 
0.2 Mpc/h. The model (smooth red solid line) seems to work 
well for small scales compared to the data (jagged lines with 
errors) . 

In the three cases, we can calculate the x^ to see how 
well the models fit to the data, as 



Figure 16. Redshift space 2-point correlation function ^((t, it) 
calculated using squares in tt — it of side=0.2Mpc. Contours are: 
0.1-50 with logarithm separation=0.6. We can see clearly the fin- 
gers of God in the line-of-sight direction 



In Fig llSI we compare the model that we obtained us- 
ing the real space correlation function just found, which in- 
cludes all the non-linear effects, once we extrapolate ^(r) 
below 0.3Mpc/h, where fiber collisions prevent us from mea- 
suring clustering. We can see that we still need to explain 
the strong elongation we see in the direction line-of-sight. 



X 



2j(C(o-, vr)reai — ^(o", ■K)ni,odei)'^/err{n, a)^ (14) 



using diagonal JK errors, which we have proved in the 
errors section that work for small a. We have seen that the 
covariance is small (less than 20%) and so we believe that 
this x^ estimation should be accurate (at least to about 20% 
accuracy) . We see how the X^ improves when we add a non- 
linear bias and variation of ai2. 

We calculate the x^ in different fitting zones, which we 
define by including in our analysis all the pixels that have 
an amplitude higher than 1.6, 2.4 or 4.4 in the third model. 
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Figure 18. Redshift space 2-point correlation function ^(ct, tt) (z=0. 15-0.47) modeled with the real-space correlation function obtained 
from dcprojoction (this one has non-linear bias dependent on scale) and an effective (Ti2=380km/s. We use Eq.JSjl to calculate the model. 
Topidata (as colors), data (as colors) -|- model (solid line), model (as colors). Bottom: the same zoomed in a 



the one that is most similar to LRG (non-linear bias and 
variation of (712). The reduced x^ (total x^ divided by the 
number of individual points) varies from 2.2 to 3.2 for the 
first model (linear bias) depending on the fitting zone; from 
1.8 to 2.8 for the second model (non-linear bias); and from 
1.2 to 1.3 for the third model (non-linear bias and variation 
in (T12). The third model represents a major improvement 
respect to the other models although the fit is not perfect 
because the reduced x^ is not equal to unity. However, the 
reduced x^ remains constant when changing the fitting zone, 
so the model is consistent at different scales. Moreover, at 
small scales, a reduced x^ of 1.2 is not rule out as can be seen 
in the first panel of Fig ll4l Apart from this effect, we think 
that the difference to a perfect reduced x^ can be attributed 
to the covariance which is small but non zero for LRG and 
has been neglected in this analysis. We should also take into 
account that in order to recover (T12 we go through many 
steps, that can increase the error at the end. As an input 
for the calculation of the model of Eq.Q we need /3 and 
the real-space correlation function £,{r), and both measures 
are estimations, they are not direct observables. We know 
that the model is perfect for dark matter, but it probably 
needs the inclusion of non-linearities to explain better the 
2-halo pairwise velocity dispersions in highly biased galaxies 
as LRG (Scoccimarro 2004). We let to future work the inclu- 
sion of more realistic, non-linear corrections, to the Kaiser 



model. As a conclusion, we see that this simple model can ex- 
plain the strong FOG without need of having large pairwise 
velocities, and the model recovers almost all the features of 
5((T, tt), a complicated mix of different effects, as explained 
in previous sections. 



5.8 Different redshift slices 

Finally, we look at the differences in the bias, defined in 
Ea. (|13|l . for different redshift slices. In Paper I we saw that 
b{z)D{z) was nearly constant as a function of redshift. This 
means that b{z) grows smoothly with redshift in agreement 
to what we see in Fig |21l The slope in the non-linear bias 
is nearly the same, so the small scale interaction between 
galaxies is nearly the same as a function of redshift, as ex- 
pected. 

We have done the same analysis for different slices and 
the models for ^(cr, n) work very well for all the cases. Plots 
are very similar to the ones in previous section. 



6 DISCUSSION 

We have studied small scales in ^(cr, tt) measured from lumi- 
nous red galaxies, which have both a strong non-linear bias 
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Figure 19. Redshift space 2-point correlation function ^(ct, tt) (z=0. 15-0.47) modeled with the real-space correlation function obtained 
from doprojcction (using Eq. JSj) and ai2 dependent on scale, as plotted in Fig llSI We include the variation of ai2 with a and assume 
that (T12 is constant along the line-of-sight for a fixed a (perpendicular distance). Top:data (as colors), data (as colors) + model (solid 
line), model (as colors). Bottom: the same zoomed in a 
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Figure 21. Comparison between the bias b{r) = \/ ^{r) / ^(r) p m 
(as in Eg. mSI l) in real space for different slices in redshift. All: 
black, z=0.15-0.3 (solid blue), z=0.3-0.4 (dashed blue), z=0.4- 
0.47 (dashed-dotted blue); z=0.15-0.34 (solid red), z=0. 34-0.47 
(dashed red) 



and are affected by large random peculiar velocities in red- 
shift space. We have compared our results for the monopole 
in redshift space ^(s), the real-space correlation function ^(r) 
and the perpendicular projected correlation function E(r) 



with Zehavi et al (2005) and Masjedi et al (2006) resuhs. 
Our analysis agree at all the scales except for very small 
scales, below IMpc/h, where we find some minor differences 
with respect to previous analysis. 

Here we also recover, for the first time, the pairwise 
velocity dispersion ai2 of LRG galaxies as a function of sep- 
aration by fitting ^(ct, tt). Our method is shown to work in 
realistic simulations. We find that the scale variation of ai2 
in LRG galaxies is similar to that in dark matter simula- 
tions. We show in section i|5.7l that a simple Kaiser model 
convolved with an exponential distribution of pairwise ve- 
locities, can explain well the complicate shape of ^{(j, n) at 
small scales, once we add the scale dependent bias and the 
scale dependent ai2- The x^ P^r degree of freedom reduces 
by over a factor of two when we allow 6 and o"i2 to change 
with scale. We notice that if we attribute all the distortion 
at small scales to the peculiar velocities, without taking into 
account the non-linear bias, velocities appear to be artifi- 
cially larger. On small scales the errors in ^(a, n) are rela- 
tively small compare to the signal, see section i]3.1l Thus the 
agreement of the ^{a, n) data to our simple modeling shown 
in Fig ll9l is quite remarkable and significant. We believe that 
this is the first time this simple ^{a, n) model is shown to 
agree in detail with small scale LRG clustering. 
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Figure 20. We plot ^(cr, tt) as in Fig |19l along the line-of-sight n for different fixed a in bins of 0.2 Mpc/h as labelled in the panels. 
Jagged solid line with errors (gray zone) is the measured LRG ^{cr, n) while the smooth line (red) is the best model (with non-linear bias 
and (Ti2 dependent on scale). Here we only show a few values of a, but the model works well in all cases. 
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